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Abstract 

We investigate D^n^ decay in perturbative QCD approach which has recently 

been applied to B meson decays. B^ -D^7r+ decay (and its charge conjugated mode) 
can be one of the hopeful modes to determine \ Vub\ since it occurs through b —>■ u transition 
only. We estimate both factorizable and non-factorizable contribution, and show that the 
non-factorizable contribution is much less than the factorizable one. Our calculation gives 
BR{W ^ D;n+) = (50 ~ 70) x 
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The determination of all the elements of Kobayashi-Maskawa matrix is important 
for the consistency check of the standard model and a search for new physics. Much 
extensive experimental efforts have been being done at B meson dedicated facilities (B 
factories) to complete the determination of KM matrix elements to third generation. \Vcb\ 
has well been determined owing to heavy quark symmetry to the accuracy of less than 
10 % error. But \Vub\ has not yet been determined so precisely|0]. The experimental error 
will be greatly reduced by the coming experiments in the very near future, while we need 
more effort to reduce the theoretical uncertainty. It is mainly due to hadronic effects 
which we do not yet have an precise method to calculate. 

So far 6 — i> M semi-leptonic decays have been mainly used for the experimental de- 
termination of \Vub\- It is interesting to investigate other modes involving non-leptonic 
decays to extract \Vub\ for the consistency checks of the experimental value of \Vub\ and 
the theoretical methods to estimate hadronic effects. More experimental information can 
be available and we can tune up the theoretical methods. We propose here 5° D~ti^ 
decay as one of good candidates to investigate. The final state D~n^ is composed of 
(sc) (ud) quark state. The b quark in meson cannot directly decay into c quark by 
W~ emission as the color quantum number is different. Also no penguin, exchange nor 
annihilation contributions exist since no qq state presents in the final state. Therefore, 
the decay occurs through b ^ u transition only, which makes B^ Dj7i~^ a good mode 
to determine \Vub\- 

B^ DjTT^ and its charge conjugate mode are hopeful from experimental point 
of view also. Both Ds and charged pion are relatively easy to be identified in the 
present experiments. Recently, BABAR and BELLE groups obtained the branching ratio, 
BR(50 -> D+TT-) = 3.2 ± 0.9 ± 1.0 x 10"^ (BABAR) g, 2.4 +J;g ± 0.7 x 10"^ (BELLE) 
1^. With increasing statistics at B factories we can expect that the branching ratio is 
fixed more precisely in the very near future. 

The B^ D^TT^ decay occurs through the effective Hamiltonian, 

Ti-eff = --^VubVcs*[ci{sb)v-A{uc)v-A + C2{ub)v-A{sc)v-A] + ih..C.) (1) 
= ^K6K/[ClOi + C202] + (h.C.), (2) 

where Ci 2 are the Wilson coefficients obtained by solving renormalization group equa- 
tions, and Oi^2 are 4-quark operators]^. We need to estimate {D'^n^\0i^2\B^) to obtain 
the branching ratio theoretically. So far, factorization ansatz^] has often been used to 
estimate this kind of 2-body decay hadron matrix elements; 

A[W^D;n^] ^ %K6K/fc2 + -^^ {D;{P2)\{sc)v-A\0){n+{P3)\{ub)v-A\B^{Pi)) 



^'uo^cs y^. , ^^^^^ 



^fDs^V^bVcs* (c2 + 4^] M|Fo(Mz),2), (3) 
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where N^ff is the effective color number and fos is the Dg meson decay constant. Pion 
mass is neglected. The B ^ n transition form factors are defined as follows; 



(7r+(P3)|^7M&|5°(Pi)) = 



(^1 + Ps), - 



'Mi; 



(4) 



where q = P\ — P3. Form factors are to be obtained from another theory or experimental 
data. The deviation of iVe// from the number of color, = 3, accounts for so-called 
non-factorizable contributions. With the amplitude given by eq.(|]) the -D^vr"*" 
branching ratio is calculated as 



BR(50 



T{B^y^ 



32n 



-GlM\ 



C2 + 



Cl 



N, 



Fo(Mb.)VI>sIK6Ks| 



(5) 



eff , 



where r = M^s/Mb- Below we first give estimations of the branching ratio based on the 
factorization ansatz by using the B ^ n transition form factor from light-cone sum rules 
and lattice QCD. Then we calculate B^ — » -D^tt^ amplitude by using perturbative QCD 
(pQCD) approach to estimate non-factorizable contribution. pQCD approach has been 
applied to estimate pion electro-magnetic form factor, B ^ D, B ^ n transition form 
factors and several decay amplitudes of 2-body decays of B meson {Dn, tttt, Ktt and K*'-f) 
The results of pQCD approach nicely agree with experimental data. The advantage of 
pQCD approach lies in the point that the non-factorizable contribution can be calculated 
based on well-established perturbative QCD technique for heavy meson decays. We show 
that the non-factorizable contribution is about 10% or less of the factorizable one in 
B DgTT, so that naive factorization estimation gives a reasonable prediction for this 
decay mode. 

In the evaluation of eq.(|]) we consider the Wilson coefficients at two scales, /x = Mb 
and Mb/2, to estimate the ambiguity coming from the choice of the scale; 





Mb 


Mb/2 


Cl(/^) 


-0.274 


-0.393 


C2(/i) 


1.12 


1.19 



The S — > TT transition form factor based on light-cone sum rules calculation can be 
parametrized as 

^o(O) 



1 



aF{q^/Ml) + hp{q^/Miy 



(6) 



with Fo(0) = 0.305, ap = 0.266 and bp = -0.752 0. With the calculation by lattice 
QCD[|TT| we make a fit for B ^ 71 transition form factor as 



FoiO) 



CFig'/Ml) 



(7) 
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Neff = 3 


Neff = 2 


Neff = 3 


Neff = 2 






Ai = MB 


= Mb/2 


/i = Mb/2 


lattice Fq 


51.1 


46.7 


53.9 


47.4 


sum rules Fq 


44.0 


40.1 


46.4 


40.8 



Table 1: BR(fiO ^ /^;7r+)/(/|,jKfcKsn 



with Fo(0) = 0.310, cp = 0.760. The estimated value of BR(50 ^ D;7c+)/{fl^^ 
IS given m Table |l|. By adopting fos = 0.241 ± 0.032 GeV obtained by taking average of 
several experimental data|T^, we can summarize the naive factorization estimation as 



BR(50 



-)/\VubV,, 




4.0 (lattice Fq) 
3.5 (sum rules Fq) 



The ambiguity from the choice of N^ff and fi is about 10% and 5%, respectively. The 
major ambiguity lies in value. 

There are factorizable and non-factorizable contribution in Djn^. When a 

gluon connects the spectator d quark and a quark in Dg meson in Fig.|l|, the contribution 
cannot be factorized as in eq.(§). The non-factorizable contribution is taken into account 



Ds 



b XO2 



71 




Figure 1: Diagrams contributing to B'^ ^ D^ti^ . Gluon is not shown. 

by changing the color number from Nq = 3 to Nf.ff in the calculation based on the fac- 
torization ansatz. However, the number to be taken as A^e// is theoretically unclear. We 
have to rely on fit of Nejf by using experimental data other than B — >■ Dsvr. But, we 
cannot simply adopt the Neff in other decays because the topology of diagrams contribut- 
ing to other decays is not necessary same as in the case of B'^ -D^vr^ decay. Here we 
calculate B'^ —>■ D^n^ amplitude in pQCD approach^ to estimate the non-factorizable 
contribution based on QCD. 

Some of representative diagrams contributing to B^ — >■ D~tt^ in pQCD approach is 
shown in Fig.Q. The point is that in 2 body decays of B meson the spectator quark has 
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Ds Ds" 




d 



(c) 

Figure 2: Some diagrams contributing to S° — > -Dj7r+ in pQCD. Gluon is shown in dashed 
hne. 




Figure 3: Momentum flow in the diagrams contributing to — > 7r+ in pQCD. 
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to obtain high momentum to form a meson with one of the emitted rapid quarks from 
b quark decay. The high momentum is carried by a gluon, so the perturbative QCD 
treatment is possible. The non-perturbative nature is put into meson wave functions. For 
more details refer the papers in IQ. There is another approach of calculating non-leptonic 
2-body decays of B meson by using perturbative QCD, the so-called QCD factorization 
methodp. One of the major differences between pQCD approach and QCD factorization 
lies in the treatment of the transverse momentum. It is taken into account in the pQCD 
approach while neglected in the QCD factorization. Without the transverse momentum 
the diagrams given in fig.^ have infrared singularities which makes the results of the 
calculation implausible. The calculation based on pQCD approach can be done in a 
similar way to calculate B Dn amplitude given in Li and Melicp[]. The assignment 
of momenta is shown in Fig.| where Pi = (Mb/V^) (1, 1, Ot), P2 = (Mb/V^) (r^, 1, Ot), 
P3 = {Mb/V2) (1 - r2,l,0T) and h = {0,Xi{MB/V2),hT), = (0, X2(MB/y2), 
ks = (xsi^l — r^){MB/ V^), 0, k^x) in light-cone coordinate. We obtain 



r[po 



G 



1287r 



B' 



■\fDsi] 



int 



M 



int 



(9) 



The factorizable contribution, fos^int, is given as 



5 int 



nl rl/A 

WttCfV^M'^ / dxidxs / hidhih2,dh2,4>B{xi^hi] 
Jo Jo 



Xas(^int) C2(tint) + 



Cl(tint) 



exp[-5'ij(tint) - S^{ti^t) 



c 



1 + 



X [{I + x-sil - r^)}Mx3) + ro{y^ _ J 
-ro(2x3 - l)4>T{x3)]h{xi,X3, bi, 63, mint) 



2X3}0p(x3) 



+ [xir^Mx3) + 2ro{l 



1 -r2 



Xi}(l)p{x3)]h{x3, Xi, 63, 61, mint) 



(10) 



where exp[— (tint) — 'S'^(^int)] is Sudakov factorg], mint = (1 ~ '^'^)Mbi ^ — ^qcd and 
tint = max(y'ximint, ^/x^rn^, l/^i, l/^s)- The parameter 6j (z = 1 ~ 3) is the conjugate 
variable of fej^. The function h{xi, X3, 61, 63, m) is defined as 



/i(a;i,a;3,6i,63,m) 



St{x^)KQ {^xix^mhi) 
X [e{hi - h)KQ (y/xsmbi) Iq {^x^mb^) 
+9{b3 - bi)Ko {y/xambs) Iq {^Jx^mbi] 



where 



SAx) 



2i+^T(c + 3/2) 

Ar(c+i) 



[x{l-x)Y (c = 0.3 -0.4), 



(11) 
(12) 



6 



which comes from threshold resummation||13|| . The non-factorizable contribution is given 

as 

A^int = -32'jTy2NcCFV^M'^ J dxidx2dx3 J bidbib2db2(pB{xi,bi)(f)Ds{x2,b2) 
X as {td) ^^^^ exp [-S {td) \b3=bi ] 



c 



xi — X2jr 



[{xi - X2){l + r )(t)A{,xz) +ro{a;3 _ ^^^^ }<pp{.X'i) 



(a;i — X2)r'^^ 
(l-r2) 



}'pT{x:i)]S[2h^P {xiM 



V f / 'y \ \ ~\ / \ f \ I 3^ \ I 2/ 2 / ^ ~1 y / \ 

+ - xi - X2 + (1 - r )x^)\(I)a{x^) - rojxs j^pia^s) 



(1 -r2) 



-ro{xi, + 



fl-r2) 



(13) 



where A^c = 3, = max(DMB, ^l/^f |Mb, y'l^ilMs, I/62), 

= XiX3(l-r2), 
-Di = {xi - X2)x^{l - r^) , 

(xi + a;2)r^ - (1 - xi - X2)a;3(l - r^) 



-^2 



The functions h^P and /i^^'' are defined as 



- 62)^^0 {DMBbi) Jo (I^Mb62) 
+9{b2 - b,)Ko {DMBb2) h (DMsb,)] 
Ko{DjMBb2) for D] > \ 
D'^AMnbo) for < J 



(14) 



(15) 



The factor S\n\ accounts for threshold resummation in non-factorizable contribution in- 



vestigated in the work by Li and UkailT^. The effects of this factor shall be discussed 
later. 

We have included twist 3 component into wave functions; 

B{P) = [/P + Mb]750b(x), (16) 
Ds{P) = ^^[^ + Mds^Ds{x), (17) 
7r+(P) = 75[/P0^(x) +mo0p(x) -mo(/' /i-t;-n)0T'(a;)], (18) 

where mo = M^/{mu + md), tq = itio/Mb, v = (1, 0, Oy) and n = (0, 1, 0^) in light-cone 
coordinate. We have adopted the following functions for B meson and pion; 
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(f)B{x,b) = Nbx (1 — x) exp 



(19) 



<t>A{x) = -^^x{l-x)[l + a2CT{l-2x) + a^cT{l-2x)l (20) 
0p(x) = -J^[l + a2,Cl'\l-2x) + a,,&J\l-2x)l (21) 
'/'t(x) = ^-^(l-2x)[l + 6a2t(10x2-10x + l)], (22) 

where Nb is normalization constant to have /q^ ^^(x, 0)(ix = /b/2-\/6, and is a 
Gegenbauer polynomiah Cl''^{x) = (l/2)(3a;2 - 1), Cl^^{x) = (l/8)(35x^ - 30x^ + 3), 
Ci/^(x) = 3x, C2^\x) = (3/2)(5x2 - 1), ^^^(x) = (15/8)(21x^ - Ux^ + 1). The param- 
eters in the pion wave functions are given in ref. fTT 



a2 = 0.44, 04 = 0.25, a2p = 30?73, = -30?73u;3, 02* = Sr^a - ^?73t^3, (23) 

with 7^3 = 0.015 and = —3.0. As for Dg meson wave function we take; 

^Ds = Nnxil -x)[l + |cf/'(l - 2x)]. (24) 

By analyzing B ^ n and B ^ D form factors we have obtained a set of parameters 
consistent with experimental data and other theoretical predictions JT3|, |TB[; 

ujb = ^'b = 0-4 for B meson and q = 0.7 for D meson. 

The Ds meson wave function is thought to be similar to D meson wave function by SU (3) 
flavor symmetry. The parameter q are varied in the numerical analysis to see their effects; 

Cd = 0.5, 0.7, 0.9. 

The numerical results on the B^ D~ti^ branching ratio is given in Table |] by taking 
c = 0.35 and 5*;^^ = 1; i-^- without the threshold resummation effect in non-factorizable 
contribution. The results show that the branching ratio of this mode is almost insensitive 
to the Ds meson wave function for a reasonable range of parameter q ||16|| . 

The effect of the threshold resummation is also investigated. The threshold resumma- 
tion parameter c in eq.(12) for the factorizable contributions is varied within the range 
consistent with the 5 — > vr form factor analysis |]T^. As for the non-factorizable con- 
tribution we take S[ll = St{xs) and S^^l = 1 following the arguments given in ref. P . 



The parameter c is taken to be the same with that for factorizable contributions for the 
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Cd 


0.5 


0.7 


0.9 




0.093 


0.088 


0.082 


BR/(/2jKfeK/) 


60 


59 


59 



Table 2: BR(i?'^ D^TT^)/\VubVcs\^ and the ratio of non-factorizable contribution to 
factorizable one for different Dg meson wave function parameter q. f^s should be given 
in GeV. 



c 


0.30 


0.35 


0.40 


(a) 


A^int/ fos^int 

BR/{fl^\VM') 


0.079 - 0.022i 


0.085 - 0.024i 


0.090 - 0.026Z 




68 


59 


52 


(b) 


■Mint/ fos^mt 

BR/{fl^\VM') 


0.069 + 0.018i 


0.074 + 0.025Z 


0.077 + 0.031Z 




66 


58 


51 



Table 3: (a) BR(i?° D^7i^)/\VubVcs\'^ and the ratio of non-factorizable contribution 
to factorizable one for different threshold resummation parameter c without threshold re- 
summation in non-factorizable contribution, (b) same as (a) with threshold resummation 
in non-factorizable contribution. 



simplicity of the calculation. Two kinds of calculations have been made with or without 
threshold resummation factor in the non-factorizable contribution. The results are shown 
in Table ^ It is found that the branching ratio varies about 15% depending on the choice 
of the parameter c. But this parameter is just a numerical convenience to fit the true 
form of the following threshold resummation factor by eq.(12); 



Sfix) 



dN J{N) 


—ioo 2m N 



l-x) 



-N 



(25) 



where a is an arbitrary real constant larger than all the real parts of the poles in the 
integrand and 

'1 /■! , 1-^^-1 /-(l--)' dX , ' 

-lK{as) 



J{N) = exp 



2 Jo 



dz- 



{i-z) A 



(26) 



with 7_ft:(as) being the anomalous dimensions [|T^, |T6[. Thus the parameter c is absent in 
the more careful (but complicated) treatment of numerical calculation where the above 
equations are used. So this ambiguity does not lead to a true theoretical error. The effect 
of the existence of the threshold resummation in non-factorizable contribution is found 
to be small in this decay. This is because non-factorizable contribution is small, 10% or 
less in amplitude, in comparison with the factorizable one in this decay mode. But the 
effects of the threshold resummation in non-factorizable part can be significant in another 
modes such as color-suppressed decay like Z}°7r°||l8 . 
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The prediction by pQCD approach is given considering the ambiguity discussed before 

as 



BR{B^^D;n+) = (50 ~ 70) X /2 jKfeK/ (27) 
= (2.4 ~ 4.6) X \VubVcs\^. (28) 

If we take the central values of parameters and the experimental data of i? ^ DsTt 
branching ratio given by BABAR and BELLE, we obtain, 

iKfcl = (3±1) X 10-^ (29) 

which is in good agreement with the value of \Vub\ obtained in 6 — m semi-leptonic 
decay[|l^]. This agreement gives a support to our treatment of B meson decays in pQCD 



approach. It also implies that the naive factorization ansatz works well also in this decay 
mode[|l^ since our calculation has shown the dominance of factorizable contribution. The 
ambiguity will be reduced when fj^g value is fixed more precisely in the future experiments. 
In the near future high statistics of B meson decay data will be available, then we can 
make pQCD prediction more precise by fitting parameters of the wave functions with rich 
experimental data. Then \Vub\ can be determined as precise as those from semi-leptonic 
decay by using the branching ratio of B^ D^n^ . 
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